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ABSTRACT 

The current state of the art in pulsar magneto spheric modeling assumes the force-free limit of 
magnetospheric plasma. This limit retains only partial information about plasma velocity and 
neglects plasma inertia and temperature. We carried out time-dependent 3D relativistic mag- 
netohydrodynamic (MHD) simulations of oblique pulsar magneto spheres that improve upon 
force-free by retaining the full plasma velocity information and capturing plasma heating of 
putative magnetospheric gamma-ray emission regions. We present the velocity structure in an 
oblique pulsar magnetosphere and the dependence of pulsar spindown power on inclination. 

Key words: MHD — stars: neutron — pulsars: general — gamma-rays: theory — methods: 
numerical — relativity 



1 INTRODUCTION 

The launch of the Fermi satellite opened a new window into study- 
ing the properties of pulsar magneto spheres, with more than 100 
confirmed gamma-ray pulsars to date (Nolan et al. 2012). Since 
gamma-ray emission can comprise a substantial (< 10%) fraction 
of total pulsar spin-down losses, gamma-ray production mechanism 
should be capable of efficiently converting magnetospheric electro- 
magnetic energy into gamma-ray radiation. 

Historically, much effort went into studying vacuum pulsar 
magnetospheres based on the analytic magnetic field solution by 
Deutsch (1955), even though it was realized early on that for pul- 
sars, which are filled with abundant plasma (Goldreich & Julian 
1969), force-free approximation is a more appropriate framework. 
Force-free approximation accounts for magnetospheric charges and 
currents but neglects plasma inertia, which is appropriate in the 
limit of high magnetization of the plasma cr = b 2 /Anpc 2 » 1, 
with b and p the proper magnetic field and density. In recent years, 
self-consistent force-free solutions of axisymmetric (Contopoulos, 
Kazanas & Fendt 1999; Gruzinov 2005; McKinney 2006b; Timo- 
khin 2006) and oblique (Spitkovsky 2006; Kalapotharakos & Con- 
topoulos 2009; Petri 2012) pulsar magnetospheres were numeri- 
cally obtained. They feature a magnetospheric current sheet that 
could be responsible for powering the observed gamma-ray emis- 
sion (Bai & Spitkovsky 2010a,b). 

Magnetohydrodynamic (MHD) approach has important ad- 
vantages relative to force-free because it includes plasma inertia, 
thermal pressure, and velocity component along the magnetic field. 
While all of these are crucial for current sheet physics and for the 
construction of gamma-ray light curves, all are neglected in force- 
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free models. MHD has no such uncertainties, paving the path for an 
improved gamma-ray light curve modeling and, ultimately, the un- 
derstanding of global structure, dissipation, and emission processes 
in pulsar magnetospheres. The first two-dimensional MHD simula- 
tion of the aligned pulsar magnetosphere was presented by Komis- 
sarov (2006). In this Letter we extend this work and construct time- 
dependent 3D relativistic MHD models of both axisymmetric and 
oblique pulsar magnetospheres and make comparisons to force- 
free models. In §2 we describe our numerical method and problem 
setup, in §3 we present our results on magnetospheric velocity, ther- 
mal pressure, dissipation, & energy losses, and in §4 we conclude. 



2 NUMERICAL METHOD AND PROBLEM SETUP 

We use general relativistic MHD code HARM (Gammie et al. 
2003 ; McKinney & Gammie 2004) including recent improvements 
(Tchekhovskoy et al. 2007, 2011; McKinney & Blandford 2009; 
McKinney et al. 2012). For direct comparison to previous studies, 
we neglect stellar gravity and carry out the simulations in flat space. 
We place the surface of the neutron star (NS) at r* = 0.2R LC , where 
^lc = c/Q is light cylinder (LC) radius, Q is pulsar angular fre- 
quency, and P = 27r/Q is the pulsar period. We use a spherical 
polar computational grid, r, 6, (p, with 6 = along the rotation axis, 
Q. We also sometimes use cylindrical radius, R = r sin 6. The grid 
extends from the NS surface, r in = r* = 0.2R LC , to r out = 200^ L c. 
The 0- and <p- grids extend from to n and to 2n, respectively. 
Grid spacing is uniform in 6- and ^-directions. The spacing is log- 
arithmic in the radial direction, Ar/r = constant, at r < r = 20/? LC , 
and becomes progressively sparse, Ar/r oc (log r) 3/4 , at r > r . 

We initialize the simulations with a dipolar magnetic field of 
a dipole moment, //, that makes an angle, or, with the rotation axis, 
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Figure 1. Meridional cut through relativistic magnetohydrodynamic magnetosphere for an aligned pulsar (model D0R2048) taken after 20 rotations; by this 
time the shown part of the magnetosphere reached steady state. Solid lines show contours of enclosed magnetic flux which represent poloidal field lines. 
[Panel (a)] Colour shows out-of-plane magnetic field component, B® (red/blue - pointing into/out of plane, see colour bar), associated with poloidal currents 
circulating in the magnetosphere. The currents flow out of the polar cap and return through the current sheet along the midplane (seen as the sharp blue-orange 
transition) and the current layer along the boundary of the closed zone (shown with thick solid line). [Panel (b)] Colour shows the logarithm of the ratio of 
enthalpy, w = p + (T + l)u g , to b 2 (red shows high and blue low values, see colour bar). This ratio is small near the star, hence the magnetosphere is highly 
magnetized (w/b 2 < 10 -3 , the magnetization is about lOx that in Komissarov 2006). However, w/b 2 exceeds unity in the equatorial current sheet, in which 
thermal pressure slows down plasma inflow into the sheet, thereby changing the velocity structure inside and near the sheet. 



Q. At the inner r-boundary (the stellar surface), we set the per- 
pendicular (to the magnetic field) 3 -velocity component, v ± , to 
enforce stellar rotation, and set the parallel 3 -velocity component 
to zero, v\\ = 0. We apply at the outer r-boundary standard out- 
flow boundary conditions (BCs), at ^-boundaries transmissive po- 
lar BCs (McKinney et al. 2012), and at (^-boundaries periodic BCs. 

In a typical pulsar the magnetization near the LC can be very 
high, cr LC ~ 10 4 . In a dipolar field, quantities drop off rapidly with 
r: p oc b oc <t oc r~ 3 . This is numerically challenging: to ensure 
force-free-like conditions at the LC, <x LC » 1, we must have a 
very high cr near the star, cr* = (Rlc/R*) 3 o~lc - 10 2 cr LC , i.e., 
much higher than cr ~ 10 2 that our code can handle at reason- 
able resolutions in 3D. For our simulations to closely resemble the 
actual pulsar magneto spheres, we ensure force-free-like conditions 
inside the LC by driving density, p, internal energy, u g , and par- 
allel 4-velocity, u\\ = yv\\, to target values, p t = b 2 IAnc 2 cr t (with 
cr t = 50 or 100), u g j = 0.2p r c 2 , and uy = 0. We do this via mod- 
ifying the variables, q = {p,u g ,u\\}, at the end of each time step 
via q = q t + (q - q t ) exp(-/cAf/r), where k - {cos 6 m , cos 6 m , 1}, 
6 m is the magnetic colatidute, At is the time step. For the driv- 
ing timescale, we choose r = at r < r\ = 0.5R LC , r = 
ToO - ri)(r 2 - r»)/[(r 2 n)(r 2 - r)] at n < r < r 2 = R L c, and 
r = oo at r > r 2 , where r = 0.005P (see also Komissarov 2006). 



3 RESULTS 

We carried out a number of simulations for different values of the 
magnetic dipole inclination angle relative to pulsar rotation axis, 
or, from 0° (aligned rotator) to 90° (orthogonal rotator). We refer 
to these as models Dxx, where xx is a measured in degrees, see 
Table 1. We indicate resolution used in a model via a suffix "R" 
followed by the number of grid cells in the r-direction; we omit 
the suffix for our default choice, = 256. Our ^-resolution is tied 
to A^ r via N e = 0.5N r , so the aspect ratio of computational cells is 
about unity, Ar : rA6 : rA(p w 1 : 1 : 2. For all MHD models we use 
an ideal gas equation of state, p g = (T - \)u g , with the poly tropic 
index, Y = 4/3, appropriate for a relativistically hot pair plasma. 
Figure 1 shows the structure of magnetic field and the ratio of 



Table 1. Model details: model name, pulsar obliquity angle (a), simulation 
resolution (given as Af r X Ng X N<p), final time of the simulation (tf, mea- 
sured in units of pulsar period, P), pulsar spindown luminosity (L, measured 
in units of aligned pulsar spindown luminosity, L a ii gn ed> ec l- 1)» f rac ti° n of 
spindown luminosity dissipated in the inner wind zone, r < 5/?lc (e). 



Name 




Resolution 


t f /P 


L/ ^aligned 


6[%] 




Relativistic MHD models (with HARM): 




D0R64 





64 x 32 x 1 


33 


0.929 


28 


D0R128 





128 x 64 x 1 


45 


0.983 


23 


DO 





256 x 128 x 1 


63 


1 


17 


D0R512 





512 x 256 x 1 


44 


1 


14 


D0R1024 





1024 x 512 x 1 


45 


0.994 


11 


D0R2048 





2048 x 1024 x 1 


22 


0.988 


8.7 


D15 


15 


256 x 128 x 128 


5.2 


1.13 


19 


D30 


30 


256 x 128 x 128 


5.1 


1.36 


17 


D45 


45 


256 x 128 x 128 


5.1 


1.64 


15 


D60R64 


60 


64 x 32 x 32 


4.8 


1.92 


27 


D60R128 


60 


128 x 64 x 64 


5.5 


1.94 


16 


D60 


60 


256 x 128 x 128 


3.3 


1.92 


13 


D60R512 


60 


512x256x256 


3.5 


1.96 


12 


D75 


75 


256 x 128 x 128 


2.9 


2.14 


12 


D90 


90 


256 x 128 x 128 


5 


2.2 


11 




Force-free models (with HARM) 






D0R64ff 





64 x 32 x 1 


7.7 


0.839 


50 


D0R128ff 





128 x 64 x 1 


30 


0.886 


47 


DOff 





256 x 128 x 1 


240 


0.914 


45 


D0R512ff 





512 x 256 x 1 


10 


0.925 


45 


D0R1024ff 





1024 x 512 x 1 


16 


0.93 


43 


D0R2048ff 





2048 x 1024 x 1 


13 


0.932 


43 


D30ff 


30 


256 x 128 x 128 


6.9 


1.27 


23 


D60ff 


60 


256 x 128 x 128 


7 


1.83 


5.8 


D90ff 


90 


256 x 128 x 128 


7 


2.11 


2 



gas enthalpy to magnetic energy in our highest resolution aligned 
model, D0R2048, which is a 2D simulation. Since our relativis- 
tic MHD models are highly magnetized, with magnetization inside 
the LC <r t = 100 » 1 (§2), they display the same generic features 
as force-free models. Unlike vacuum magnetospheres, which lack 
charges and currents, both relativistic MHD and force-free mag- 
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Figure 2. Radial profiles of angle-integrated Poynting flux. The Poynting 
flux is normalized to the analytic approximation for the spin-down of an 
aligned rotator, see eq. (1). The solid lines show total energy flux and dashed 
lines show Poynting flux. For aligned pulsar, electromagnetic fluxes at dif- 
ferent resolutions are shown. From bottom to top, the figure shows models 
D0R64-D0R2048 (red), D30 (green), D60 (blue), D90 (magenta lines). 



netospheres readily supply charges and currents necessary to sus- 
tain corotation of the magneto sphere with the pulsar. These currents 
flow out of the polar cap of the pulsar and return in the equatorial 
current sheet, which appears in Fig. 1(a) as a sharp orange-blue 
transition as changes sign across the sheet. The currents in the 
magneto sphere and in the sheet cause the magnetosphere to open 
up and form a radial Poynting-flux-dominated wind (Michel 1973; 
Contopoulos et al. 1999; Gruzinov 2005, 2008; McKinney 2006b; 
Spitkovsky 2006; Kalapotharakos & Contopoulos 2009; Li et al. 
2012b; Parfrey et al. 2012). However, not all of this wind reaches 
infinity: part of it enters the current sheet and heats it, possibly caus- 
ing the sheet to produce high-energy emission (Bai & Spitkovsky 
2010a,b). Whereas force-free approximation neglects thermal pres- 
sure, Figure 1(b) demonstrates that in MHD, reconnection heats up 
the current sheets so they are dominated by thermal pressure. 

We confirmed that our results are numerically converged. For 
this, we carried out simulations of aligned pulsar magnetospheres 
over a wide range of resolution (the first 6 models in Tab. 1). Red 
lines in Figure 2 show radial profiles of energy flux for some of our 
axisymmetric models (see the legend at the bottom right of Fig. 2). 
The total energy flux is independent of resolution and distance and 
agrees to within 1% with the spindown losses found for aligned 
pulsars (Gruzinov 2005; Spitkovsky 2006; McKinney 2006b), 

n 2 & 



-^aligned — 



(1) 



We quantify the amount of dissipation in the wind zone as 
a fraction of total energy flux dissipated in the interval r < 5R LC 
via 6 = 1 - L EM (5R LC )/L(R LC ). Table 1 and Fig. 2 show that e 
monotonically decreases with increasing resolution, e oc Ny ll?> = 
(2N e )~ 1/3 : our axisymmetric MHD magnetospheres asymptotically 
(in the limit of infinite resolution) become dissipationless. This is to 
be expected: the level of current sheet dissipation is controlled by 
the magnetospheric resistivity, which in our approach is determined 
by the numerical resolution (see also Lyutikov & McKinney 201 1). 

Our relativistic MHD magnetosphere is similar to the one 
obtained using the low-dissipation force-free code by Spitkovsky 
(2006): most of the field lines that cross the surface of LC open 



up to infinity, with only a small fraction of them entering the cur- 
rent sheet, where they dissipate a vanishingly small fraction of pul- 
sar spindown energy. Note that due to high numerical dissipation, 
standard force-free codes often reach a very different solution: most 
poloidal magnetic field lines close through the midplane, where 
they dissipate most of pulsar spin-down energy. Komissarov (2006) 
noted this effect, and we do as well when we use the force-free ver- 
sion of HARM (McKinney 2006a; Lyutikov & McKinney 2011). 
Namely, we find very high levels of dissipation, e ~ 50%, that 
do not decrease with increasing resolution (see Tab. 1). Such ab- 
sence of decrease is unphysical: had the models included thermal 
pressure produced by current sheet dissipation, that pressure would 
have slowed down the inflow into the current sheet, suppressed 
the reconnection, and led to lower values of e (Komissarov 2006). 
Recently, Gruzinov (2008, 2011a,b,c, 2012) considered resistive 
force-free simulations of pulsar magnetospheres and argued that 
magnetospheric current sheet dissipates up to 50% of pulsar spin- 
down energy and that an increase in resolution does not change 
the dissipation rate. Our results suggest that the results of Gruzinov 
(2008, 201 la,b,c, 2012) for axisymmetric pulsars are dominated by 
the large numerical resistivity of their numerical scheme and con- 
verge to the unphysical, dissipative force-free solution. 

Why does the force-free version of HARM (McKinney 2006a; 
Lyutikov & McKinney 201 1) and many other force-free codes (e.g., 
Komissarov 2006; Gruzinov 2011c, 2012; Petri 2012) show such 
high levels of dissipation? To handle discontinuities in the flow, 
force-free HARM uses a Lax-Friedrichs Riemann solver that is not 
specialized to treating current sheets, and hence spreads the sheet 
over several grid cells. However, no force-balance inside the sheet 
is possible since in force-free there is no thermal pressure to slow 
down reconnecting fields. Unless one prescribes a zero velocity of 
inflow into the sheet (McKinney 2006b), the lack of force-balance 
across the sheet in force-free causes rapid reconnection (Komis- 
sarov 2006). In contrast, the force-free scheme by Spitkovsky 
(2006) can treat current sheets as true unresolved discontinuities, 
so reconnection in the sheet is minimal. As evidenced by our MHD 
results, we believe that in the limit of low reconnection the low dis- 
sipation force-free solutions as in Spitkovsky (2006) are more rep- 
resentative of the physical pulsar magnetospheric shape than dissi- 
pative force-free solutions with uncontrolled numerical reconnec- 
tion rate (e.g., Gruzinov 2011c, 2012). 

We now consider oblique models, applicable to the majority 
of pulsars. We present results for magnetization cr t = 50 (results 
at cr t = 100 are similar and are not shown). Figure 3(a) shows the 
ft-Q plane for our oblique model D60, with a = 60° obliquity: 
electromagnetic quantities in our relativistic MHD models repro- 
duce, as expected, major features of oblique force-free solutions 
of pulsar magnetospheres (see Fig. 2(a) in Spitkovsky 2006 and 
Fig. 1(a) in Li et al. 2012b), such as the formation of closed and 
open field line zones and the undulating equatorial current sheet. 
Please see Supporting Information for movies. Additionally, MHD 
models provide crucial information about plasma properties, e.g., 
velocity and temperature in the current sheet, which are needed for 
light curve computation but are missing from a force-free descrip- 
tion. Although resistivity in our scheme is numerical, the reconnec- 
tion displays physical characteristics. Figure 3(b) shows that ther- 
mal and magnetic pressures become of the same order inside the 
current sheet, just as in our axisymmetric models (see Fig. 1(b)). 
Therefore, the thermal pressure is dynamically important in the cur- 
rent sheet and does affect the current sheet structure, e.g., velocity. 

Figure 3(c) shows the radial dependence of various 4- velocity 
components along a line y = z = that intersects the current sheet, 
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Figure 3. Slices through a relativistic magnetohydrodynamic simulation of an oblique pulsar magneto sphere (a = 60°, model D60) taken after 3 rotations. See 
Supporting Information for movies. Solid lines show field lines as traced in the image plane. Panels (a) and (b) show slices in jl-Q, or x - z, plane. [Panel (a)] 
Shows out-of-plane magnetic field component, B® with colour (red - into plane, blue - out of plane, see colour bar). Similar to the aligned pulsar (Figure 1), 
currents flow out of the polar cap and return through the current sheet. [Panel (b)] Shows quantity log 10 (w/& 2 ) in color. It is low near the star, indicating a 
highly magnetized flow (see the colour bar), and high in the current sheet, indicating the importance of thermal pressure support. [Panel (c)] Runs of 4-velocity 
components vs. radius along the y = z = line for the simulation, magnitude of 4-velocity, u = yv (green solid), 4-velocity component parallel to the magnetic 
field, u\\ (dash-dotted cyan line), and the analytic split-monopolar force-free solution for the magnitude of 4-velocity, w ff ' mono = QR (dashed magenta line, 
see legend). The split-monopolar solution systematically under-estimates u and neglects u\\. The location and thickness of the current sheet is indicated by the 
yellow stripe. [Panels (d)-(f)] Show colour maps of u in three orthogonal slices (see legends). In the sheet, which is the putative source of gamma-ray photons, 
the 4-velocity (and its individual components, not shown) clearly undergoes rapid changes, which affect directionality and beaming of emergent radiation. 



indicated by a vertical yellow stripe at x « 2.2R LC . Clearly, the 4- 
velocity component parallel to the magnetic field in the bulk of the 
flow, except near the sheet, is positive, u\\ > 0, hence, plasma slides 
along magnetic field lines predominantly away from the star. Inside 
the current sheet, all velocity components undergo a smooth varia- 
tion that manifests itself as a jump across the sheet in all velocity 
components, e.g., the magnitude of 4-velocity, u = yv, undergoes a 
30% drop across the sheet, due to reconnection-induced inflow of 
fields and plasma into the sheet. This behavior is not included in 
the split-monopole model of an oblique rotator (Bogovalov 1999), 
in which the magnitude of 4-velocity increases linearly with cylin- 
drical radius, w ff mono = QR (e.g., Narayan et al. 2007). Note that 
inside the current sheet, the velocity components in the simulation 
appear to pass through the split-monopole values (u = £IR, u\\ = 0), 
which suggests that in a volume- averaged sense the split-monopole 
model might give a reasonable description of velocity in the cur- 
rent sheet in the presence of reconnection. The lower row of panels 
in Fig. 3 shows three orthogonal cross-sections through our fidu- 
cial oblique model, D60. Clearly, the jump in u across the sheet is 
a universal property and indicates an inflow into the current sheet 
from both sides due to the magnetic reconnection. 

Figure 2 shows the radial profiles of total and Poynting (elec- 
tromagnetic) energy fluxes for models at different values of pulsar 
obliquity. The total energy flux, L(r), is conserved to better than 
3% accuracy in all cases. We confirmed that our oblique models 



are numerically converged: a twofold increase in resolution affects 
the spindown rate by less than 2% (compare models D60R128, D60 
and D60R512 in Tab. 1). Based on this, we conservatively estimate 
that the accuracy of our spin down energy loss measurements is 
^5%. Interestingly, Table 1 shows that changing resolution for our 
tilted model D60 leads to much smaller changes in the amount of 
magneto spheric dissipation, e, than in our axisymmetric models. 
This suggests that the 3D motion of the current sheet through the 
numerical grid increases the dissipation in our scheme. 

Table 1 shows that the amount of magnetospheric dissipa- 
tion in our force-free HARM models, D0ff-D90ff, dramatically de- 
creases as inclination angle increases. Namely, dissipation is un- 
acceptably high, e > 25%, at low inclination angles, a < 30°. 
However, the dissipation becomes much smaller, e < 5%, at higher 
inclination angles, a > 60°. This low level of magnetospheric dis- 
sipation (in fact, even smaller than in our MHD models) for highly 
inclined force-free HARM pulsars indicates that force-free HARM 
models D60ff and D90ff provide a good description of the elec- 
tromagnetic part of pulsar magneto sphere. Our MHD models do 
not show such strong trends of e vs. a and are applicable at all 
inclination angles. It is likely that the improved dissipation proper- 
ties of force-free schemes at large inclinations have to do with the 
larger fraction of the current in the sheet that is carried by the dis- 
placement current for higher obliquity. Force-free schemes without 
conduction currents become vacuum-like in the sheet region, and 
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Figure 4. Pulsar energy losses vs. inclination angle, a, for our MHD 
models, D0-D90 (red squares), and force-free HARM models, D0ff-D90ff 
(black circles). The losses are normalized by the analytic approximation for 
the aligned rotator, ^aligned (see eq. (1)), and are in excellent agreement with 
Spitkovsky (2006). 



and gamma-ray light curves (Li et al. 2012a,b; Kalapotharakos 
et al. 2012a,b). This will also allow studies of plasma accumula- 
tion and plasmoid formation near the Y-point that can explain pul- 
sar glitches and associated changes in pulsar braking indices (Con- 
topoulos 2005; Bucciantini et al. 2006; Spitkovsky 2006). 
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this may suppress reconnection there. This is not the case for MHD 
schemes which still have to include the plasma in the current sheet. 

Pulsars at high obliquities lose larger amounts of energy than 
at low obliquities (Fig. 4). Both in our MHD and force-free HARM 
models the spin-down power, L, is well-described by the analytic 
fitting formula, L/L aligned = k\{\ + k 2 sin 2 a), due to Spitkovsky 
(2006), albeit with slightly different values of numerical factors, 
k\ - \ and k 2 - 1.2. Our force-free HARM spin-down values 
are within 5% of the MHD values (Tab. 1), indicating an excellent 
agreement between the two approximations at high obliquities. 



4 CONCLUSIONS 

We obtained axisymmetric and oblique pulsar magnetosphere solu- 
tions using time-dependent relativistic MHD equations in 3D. We 
used a conservative relativistic MHD formulation that allowed us to 
account for resistive heating and thermal pressure support in mag- 
netospheric current sheets, both of which are important for obtain- 
ing numerically converged solutions (see § 3). Our solutions are 
highly magnetized, with cr = b 2 /AjTpc 1 50-100 near the LC, and 
are, therefore, essentially force-free. We verified that the electro- 
magnetic spin down power in our relativistic MHD models quan- 
titatively agrees with force-free models. Our MHD models gener- 
alize force-free solutions by providing crucial information about 
fluid motions that is missing from a force-free description: plasma 
density, pressure, and the velocity component parallel to the di- 
rection of magnetic field, u\\. Knowing this information is required 
for computing the beaming and phase of gamma-ray emission that 
emerges from the current sheet (Bai & Spitkovsky 2010a,b). These 
calculations will be presented in an upcoming publication. We note 
that while MHD models provide a complete description of plasma 
motion along the open field lines, we still have to switch to a force- 
free-like description inside the LC to avoid mass and internal en- 
ergy build-up on the closed field lines (see § 2). 

Our work considered a non-relativistic outflow from the sur- 
face of the NS (u\\ = 0). We plan to investigate if an ultra-relativistic 
outflow from the surface (u\\ » 1) can cause ultra-relativistic veloc- 
ity inside the current sheet on the scales of LC. As magnetospheric 
conductivity can vary, possibly influenced by the amount of mag- 
netospheric plasma supply (Li et al. 2012a,b), resistive relativistic 
MHD codes should be developed (Komissarov 2007; Palenzuela 
et al. 2009; Dionysopoulou et al. 2012) and used to study physical 
resistivity effects on the structure of magnetospheric current sheets 



SUPPORTING INFORMATION 

Additional Supporting Information is available for this article. 
Movie files. Movies of model D60: movie 1 (link) and 2 (link). 
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